Power laws in a 2-leg ladder of interacting spinless fermions. 
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We use the Density-Matrix Renormalization Group to study the single-particle and two-particle 
correlation functions of spinless fermions in the ground state of a quarter-filled ladder. This ladder 
consists of two chains having an in-chain extended Coulomb interaction reaching to third neighbor 
and coupled by inter-chain hopping. Within our short numerical coherence lengths, typically reach- 
ing ten to twenty sites, we find a strong renormalization of the interchain hopping and the existence 
of a dimensional crossover at smaller interactions. We also find power exponents for single-particle 
hopping and interchain polarization consistent with the single chain. The total charge correlation 
function has a larger power exponent and shows signs of a crossover from incoherent fermion hopping 
to coherent particle-hole pair motion between chains. There are no significant excitation energies. 

PACS numbers: 71.10.Pm, 71.27.+a, 05.10.Cc 
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I. INTRODUCTION 



The nrftlPfffi^nnn °^ quasi-one-dimensional 
conductorsaolmBijOoQ has shown that there are 
dimension-specific aspects not observed in conventional 
three-dimensional solids. Aside from the interaction 
dependent power law behavior of single-particle and pair 
response functions and the well documented spin-charge 
separation, there is the renormalization of the transverse 
hopping whose impact on the description of real materi- 
als is much debateaOO. As discussed in length in Ref. 
|3|, strong Coulomb interactions can dramatically reduce 
the effective value of the transverse hopping and retard 
the dimensionality crossover from a one-dimensional 
(ID) to a two- or three-dimensional conductor. The 
simplest theoretical testing ground for this idea is a 
2-leg ladder consisting of interacting spinless fermions on 
two chains coupled by a transverse hopping t±. It is in 
principle possible to study the putative renormalization 
of t±. This has been done using pSfarious approaches, 
among which exact diagonalizationq. momentum-space 
renormalizationiij, and bosonizationOOEjO. In all 
these papers, renormalization of the interchain hopping 
is confirmed. What we propose is a numerical calculation 
of this two-chain problem at quarter-filling 
efficient density-matrix renormalization grouj 
(DMRG) in order to directly measure power law expo- 
nents and the effective value of the interchain hopping. 
Although the DMRG has recently been triedca on 
spinless fermions, the analysis focused on a half-filled 
ladder and the nearest neighbor current correlations. 



We shall first present the model Hamiltonian we shall 
be using throughout and, second, the proposed DMRG 
procedure. Thirdly, we shall validate our approach on the 
single chain situation. We thereafter present the results 
for two chains and discuss the results in the light of the 
various theoretical treatments. A brief summary follows. 



II. HAMILTONIAN 

We shall ,wse the model Hamiltonian proposed by Cap- 
poni et aZll3 for two quarter-filled chains of spinless 
fermions that interact within each chain through a finite 
extent Coulomb potential and can hop between chains 
through the hopping term t±. The Hamiltonian is 



H 



fj + h.c)j 
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V(r) n j+r</ 3 rij,p 
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where Cj,p annihilates a fermion at site j (j = 1, ...,N) 
on chain /3 (/3 = 1, 2), rij^ is the occupancy at the same 
site, and V(r) — 2V/ (r + 1) is the intra-chain interaction 
between first, second, and third neighboring sites (r = 
1,2,3) with V as the interaction strength. We have set 
the intra-chain hopping element equal to one. 

We have chosen an interaction to third nearest neigh- 
bor because the work of Capponi showed that the single- 
fcrmion exponent a, characterizing the long-range single- 
chain inter-site transfer function 



Cl(j,r) 



) oc r 



,-(!+<*) 



(2) 



can become very large (a < 1.5 for V < 6). This power 
exponent is responsible for the perhaps better known sin- 
gularity in the momentum distribution at the Fermi level 
of Luttinger liquids. In the limit of small a, one has 
[n(k) — n(fci?)] ~ \k — kp\ a sign(kp — k). Large values 
of a will be easily observed and are expected to lead to 
much more important effects on the effective value of t±_ . 
Large values of a are also synonymous with strong varia- 
tions in the stiffness K. The two are related through the 
relation 

1 



o 



(K + l/K-2) 



(3) 



for spinless fermions on a chain. Consequently, the power 
law exponents of the various response functions, which 
are related to K, will also be strongly affected. 
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III. DENSITY-MATRIX RENORMALIZATION 
GROUP 

The exact diagonalization of Eq. (|l|) by CapponiEl was 
for short chains of up to twenty sites. Needless to say that 
some sort of extrapolation procedure, finite-size scaling 
in this case, was needed to obtain ground state infor- 
mation in the thermodynamic limit. We have chosen to 
use the DMRG since much longer chains can be studied. 
This, in principle at least, should take the system much 
closer to the thermodynamic limit and do away with the 
requirement of performing a finite-size scaling analysis. 

Another shortcoming of short chain lengths has to 
do with, . a ij 'cl imensionality" crossover in the interchain 
hoppingmL-ffij. For temperatures or frequencies larger 
than approximately \t±\, the chains do not "see" the in- 
terchain hopping, which is incoherent or diffusive, and 
they are approximately independent. But in the oppo- 
site situation, the chains are tightly coupled and they 
form bands having transverse dispersion. Let us illus- 
trate this in the situation of quarter filling for V = 
and an even number of fermions. An exact solution to 
two coupled chains is available. The states are labelled 
by k m = nm/(N + 1) where 1 < m < N and have en- 
ergy E±(k, n ) — — 2cos(fc m ) ± \t±\. For t± = 0, all levels 
up to m = nip = N/4 are filled with N/2 fermions. 
There is no interchain hopping. As increases, this 
remains so until E + {k mF ) — E_(k mF+ i), that is until 
\t±\ w (irvp /2)/(N + 1). Here vp is the Fermi velocity 
equal to \pl in our units. At this point there is a sudden 
change in interchain hopping since the two top levels be- 
low the Fermi level are E_ states. The total interchain 
hopping energy E± is now —2 \t± \ . The next jump oc- 
curs at \t±\ fa 3{nv F /2)/(N + 1) when E + (k mF -i) = 
E-(k, nF+ 2), after which E± = — 4|ij_|. At a given 
\t±\ , the jumps occur at N [(ttvf /2t±)(2p — 1) — 1] 
for p = 1,2..., when E± = -2p\t ± \ . Taking = 0.1 
for example, (nvF/2t±) = lOn/^/2 ~ 22. This is a 
large value. It is therefore difficult to attempt finite- 
size scaling or the DMRG under such conditions. One 
can only hope of reaching the thermodynamic limit for 
N 3> (ttvf /2t±). This behavior is surely attenuated in 
the presence of the Coulomb interaction which scrambles 
the spectrum. But short chains remain unpredictable be- 
cause of the discrete energy spectrum. Thus the longer 
chain lengths obtainable with the DMRG would circum- 
vent this potential numerical distortion. In the event 
that t± renormalizes to much smaller values than the 
bare one, this crossover phenomenon might even prove 
cumbersome to the DMRG. In order to avoid potential 
problem we chose-|to,|Use the finite system algorithm pro- 
posed by WhiteE3c30, projecting out the ground state 
of the superblock. At a given V, we started with the pro- 
cedure with largest value of we wished to consider, 
0.5 in all cases, and then gradually decreased its value 
using the previous solution as a seed. For each set of pa- 
rameters, the iterations stopped when the discontinuity 
in the ground state energy and the excitation energy E x 



at mid-course, when all block information has just been 
refreshed, were judged acceptable. This was typically for 
3 iterations. We used open boundary conditions since 
periodic boundary conditions lead to unacceptably large 
truncation errors. 

Let us lastly comment on the number of central sites to 
use in the DMRG algorithm. The long-range character 
of the Coulomb interaction complicates the calculations. 
For two inner double-sites (two sites on each chain), the 
computation resources (execution time and memory re- 
quirements) scale roughly as 4 2 (Nb) 4 where Nb is the 
number of sites in each of the side blocks. This comes 
from the coupling of the two blocks through V(r = 3). 
If one instead chooses to have three inner double-sites, 
the blocks no longer couple and the resources scale as 
4 3 48(A B ) 3 which comes from the coupling of the inner 
sites to each block. We have found that the resources 
are similar for Nb ~ 100 in qualitative agreement with 
this crude analysis. We have used two inner sites for the 
values of Nb = 42, 64, 96 and three for Nb = 128. 

IV. SINGLE CHAINS 

It is of utmost importance to test our DMRG proce- 
dure on simpler single-chain problems. There are two 
delicate aspects that need to be validated, both linked 
to the open-ended boundary conditions. The first one 
has to do with the value of N that can be chosen for a 
specific band filling. The second one concerns the numer- 
ical treatment that must be done on the data in order to 
generate information for infinite-length chains. 

A. Chain length and band filling 

The sensitivity to open ended boundary conditions 
can be illustrated for a chain of spinless fermions with 
an interaction extending only to nearest neighbor sites 
(V(r) — V5 r ,i) near half filling. The ground state and 
the excitation energy are completely different for the 
N = 2N f and N = 2N f - 1 situations, where N f is the 
number of fermions. Fig. |l| shows the ground state exci- 
tation energy in both situations for a calculation with 
Nb = 42 block states and N = 151,152. Although 
Nb seems small, the truncation error was nevertheless 
smaller than 3 x 10~ 7 for an open ended chain. A gap 
develops for the case N = 151 and Nf = 76 but not for 
the other. As for the ground states, they show a site 
occupancy rij — n + A m cos(irj — 6j) that is alternating 
between a large and a small value. This basic pattern is 
to be expected for a broken symmetry state with a re- 
pulsive interaction. But while the modulation phase 9j 
is a constant for N = 151, it varies for N = 152. We 
find 6 ~ irj/N in a calculation where the interaction is 
introduced right away in the first iteration of the finite 
size algorithm (sudden turn on) but 9 ~ 3nj /N when a 
first set of iterations is done with V = and then V is 
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FIG. 1: Ground state excitation energy of a half-filled chain of 
spinless fermions with nearest-neighbor interaction strength 
V. The number of sites is 151 (circles) and 152 (triangles). 



introduced in the second set of iterations, using the first 
one as seed (gradual turn on) . The ground state energy is 
lowest in the latter situation. The alternating occupancy 
in the ground state can be understood by looking at the 
large V limit when one might expect the fermions to seg- 
regate on alternating sites | ...1010...). It is the boundary 
conditions that will determine the modulation phase. For 
N = 151, one would expect 1 101.. .101) to be the stablest 
situation, with 6j — 0, and this would explain a uniform 
modulation and the excitation gap. But for TV = 152, 
the chain will spontaneously create a kink soliton (phase 
shift of 7T from end to end) |101...010) =*> |101...001) 
which can then redistribute itself and lead to a gapless 
excitation situation. This is found for the sudden turn 
on. However, it is also possible for the chain to generate 
additional kink-antikink excitations. The results for the 
gradual turn on confirm this and show that this situation 
is more stable. Just how many kink-antikink pairs would 
be generated is impossible to figure out with the DMRG. 
Curiously, the DMRG seems unable to yield an unam- 
biguous ground state when N — 2Nf for open boundary 
conditions. 

The single-chain Hamiltonian we have just been study- 
ing is akin to the XXZ problem for a spin 1/2 chairCJ. 
This spin Hamiltonian can be transformed, using the 
Wigner-Jordan transformation, into our problem with 
J z = V and J x = J y = 2 except for end terms 
^V(ni + n^v), involving the first and last sites, that oc- 
cur for open boundary conditions. The XXZ chain is 
known to develop a gap for J z > J x , that is for V > 2. 
For N = 2Nf, the ground state is degenerate, 1 101...010) 
and |010...101) , has no soliton because of the end site 
repulsion and has a gap. It is the pressure applied to 
the fermions by the end terms that insures the presence 
of an excitation energy and a uniform modulation am- 
plitude in this situation. In the case of our ladder with 
N = 2Nf — 1, it is the shorter length that adds extra 
pressure to the fermions and similarly leads to a gapped 



FIG. 2: Ground state excitation energy of a quarter-filled 
chain of spinless fermions with third-neighbor interaction 
strength V — 6 as a function of the reciprocal of the chain 
length. 

situation. 

In view of this dichotomy with respect to occupation 
near half filling, it is legitimate to ask if this sensitivity 
persists near quarter filling. To this end, we did a lim- 
ited incursion with gradual turn on at V = 6, Nf = 38, 
N = 149, 150, 152, and for first neighbor (r = 1), second 
neighbor (r = 1,2), and third neighbor (r = 1, 2, 3) inter- 
actions. The ground state energy per fermion for a given 
interaction range decreases slightly with N. The effective 
constraining "pressure" when N — &Nf — m (m = 1, 2, 3) 
can explain this. The ground state excitation energy E x 
remains small going from 0.046 to 0.11 as the range in- 
creases and is insensitive to N. Judging from Fig. [l], this 
is not a significant gap and is due to the finite length 
of the chains. Fig. || shows the excitation energy as a 
function of the inverse of the chain length for a third- 
neighbor interaction strength V = 6 and Nb = 42. The 
extrapolated gap for N — > oo is indeed negligibly small. 
There is, however, a variation of the modulation of the 
site occupancy of the form rij = n + A m cos(ttj/2 — 9j). 
Indeed, for N = 152, we find 6j = (3nj)/(2N). This 
situation is the generalization of the one seen above for 
half filling and gradual turn on condition, the quantum 
of phase shift being 7r/2 instead of tt. What this 8j means 
is that the Fermi momentum is downward shifted from 
its exact quarter-filled value kp = (tt/4)(1 — 3/N) due to 
fermions being pushed to the ends. The question spon- 
taneously arises as to any possible detrimental effect of 
such modulation on the correlation functions. We shall 
answer this in the following subsection. 



B. Data processing 

We have just found that a modulation in the site occu- 
pancy results from the open boundary conditions which 
pin fermions at the end site and lead to a broken sym- 
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FIG. 3: Doubly averaged single-fermion transfer function for 
a quarter-filled chain normalized to its value at the maximum 
as a function of the site of origin j. The plots stop at j — 
(N - r)/2. We have chosen N = 152, N B = 42, V = 6, and 
three different transfer distances (see legend). 



metry state. This modulation obviously makes it non 
trivial to get occupancies (the n) or correlation func- 
tions resembling those for periodic boundary conditions 
or an infinite chain. Taking the occupancy as an exam- 
ple, for the quarter-filled situation, it can be seen that 
the modulated part A m cos(7rj/2 — 9j) can be made less 
annoying by averaging over the natural four site cycle. 
Indeed, if we define for instance Hj — 4 2j=o n i+i = 
n - A m [cos (wj/2 - Qj) + sin(7rj/2 - 0j)} (3n/4N), one 
immediately sees that the modulation is reduced by 
a factor of order N~ 1 . A second averaging, fij = 

(V 16 )X)i=o£ro=o%-M+m> would further reduce it by 
another iV _1 . One can thus, for all practical purposes, 
remove the effect of the modulation for long chains, 
thus our choice of long chain lengths. We have used 
throughout what follows double averaging with interac- 
tions reaching to third neighbor at quarter filling. 

But, unfortunately, the open boundary conditions pro- 
duce yet another deformation. The smoothed quantities, 
like n, are not global quantities but rather local ones. 
They vary along the chain, the more so the closer a site 
is to the ends. We can illustrate this by lo oking at the 
profile of the single-fermion transfer function CT(j, r) de- 
fined in Eq. (||) using double averaging, N — 152, V = 6 
and N B = 42. The plots are for 1 < j < (N - r)/2 
since the results are symmetrical about this last value. 
Fi g. | ^ shows this fu nction for r = 1,19,54, normalized 
to C1((N — r)/2,r). The healing distance increases with 
r and is at the scale of r. It is quite obvious that the 
ends can have dramatic effects at the larger values of r. 
We have used this optimal positioning at j = (N — r)/2 
in all calculations reported below of quantities dependent 
on r. 

But are averaging and optimal positioning sufficient 
to obtain results corresponding to the thermodynamic 
limit? We can check this by extracting the power law 



exponent of the inter-site transfer function Cl(r) = 
C1((N — r)/2,r). In order to do so, one needs to have 
an analytic function with which to fit the data for CT(r). 
From the known form of the single-particle Green's func- 
tion of a Luttinger liquidtm, this function would read as 



Cl(r) 



nsin(/c^ r) 
k F r 



1 + (r/A) 1 



-i/2 



(4) 



where k F is the Fermi momentum. This form was ob- 
tained from bosonization and A is a characteristic cutoff 
parameter. This form, however, has to be modified to 
account for the DMRG procedure on a finite lattice. We 
find that the gap associated with a finite N introduces a 
numerical coherence length £ in much the same way that 
a temperatur|6-|T cc would. We thus propose the 

substitutiorffiJ r — > £sinh(r/£) and thus the following 
form 



Cl(r) = 



1 



n sin(fcF r) 
k F £smh(r/£) J \l + (£sinh(r/<e)/A)< 



a/2 



(5) 



We calculated Cl(r) for a chain having N = 152, V = 
1, 2. ..6, and Nb — 42. The truncation error was less than 
1CP 8 . The fitted parameters, in the range 1 < r < 100, 
are in Table |. It is seen that /cp ~ (7r/4)(l — 3/N) as 
expected. The passer law exponents a are those found 
by Capponi et al£3 although they seem systematically 
larger by 5-8 %. The coherence lengths are satisfacto- 
rily quite large, giving added credibility to [Opr fitting 
function. The Pearson correlation coefficientcj between 
£ and E~ Y is 0.986 indicating that the two parameters 
are strongly correlated. Finally, the cutoff A is of or- 
der one, the lattice parameter, as one would expect for 
fcrmions on a lattice. It can thus be concluded that the 
fitting function, Eq. (||), is quite satisfactory, with an er- 
ror margin of order 5 %, considering the large number of 
adjustable parameters. The thermodynamic limit is thus 
recovered albeit slightly handicapped by a numerical co- 
herence length £. Now how does the 2kp charge density 
fluctuation response function turn out? We also calcu- 
lated the correlation function C2(r) = C2((N — r)/2, r) 
where 



C2(j,r) = ((nj +r - (n j+r )) (n 3 - (rij))) 



(6) 



This function measures the correlation between occu- 
pancy (proportional to charge) fluctuations on site j and 
(j + r). It has the advantage of getting at the true fluc- 
tuation correlations in a broken symmetry state. There 
are, however, many wavemunber contributions to two- 
fermion Green's functionstM. One expects q = and 
4fcp contributions aside from the sought 2k F correla- 
tions. We observe that C2(r) has a fast oscillating part 
Bf(r) and a slow modulation amplitude A s (r), such that 
C2(r) ~ A s (r)Bf(r). What we did was to exponenti- 
ate the data C2(r) exp (— In (|A s (r)Q), do a fast Fourier 
transform, remove the unwanted contributions, and un- 
exponentiate back the remaining 2k F contribution. For 
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TABLE I: Various parameters calculated for a quarter-filled TABLE II: Various parameters calculated for a quarter-filled 
chain at different values of V and for Nb = 42.. chain at different values of V and for Nb = 10.. 



V 


E x 


k f 


a 


K 


i 


A 


K c 


d c 


V 


E x 


a. 


t 


K c 




d c 


1 


0.049 


0.78 


0.106 


0.63 


81 


1.5 


0.67 


100 


1 


0.20 


0.10 


19 


0.65 


30 


22 


2 


0.064 


0.78 


0.27 


0.49 


76 


1.2 


0.52 


87 


2 


0.23 


0.26 


15 


0.50 


31 


19 


3 


0.076 


0.78 


0.49 


0.39 


65 


1.2 


0.42 


69 


3 


0.24 


0.48 


16 


0.41 


35 


22 


4 


0.088 


0.78 


0.76 


0.31 


59 


1.3 


0.35 


59 


4 


0.26 


0.75 


17 


0.33 


34 


21 


5 


0.101 


0.77 


1.08 


0.25 


53 


1.4 


0.29 


53 


5 


0.27 


1.10 


17 


0.26 


32 


21 


6 


0.114 


0.77 


1.50 


0.21 


47 


1.5 


0.20 


43 


6 


0.29 


1.50 


15 


0.21 


32 


21 



spinlcss fcrmions, the charge correlation function at 2kp 
behaves like r~ 2K where K is the stiffness defined in Eq. 
(||). The 2kp filtered data could best be fit by the ana- 
lytical form 



C2(r) = Ccos(2fc F r + v3 c )r" 2Kc exp(-r/4) 



(7) 



which has a coherence length d c that is related to the 
excitation gap in a finite chain. The resulting power law 
exponent and coherence length are also shown in Table Q. 
K c is within 5 % of the calculated values of K obtained 
by inverting Eq. (||) . d c seems closely correlated to £ and 
both to E~ x . The Pearson correlation coefficient between 
d c and £ is 0.995 and 0.984 between d c and E~ l . The 
two coherence lengths £ and d c enter the fitting functions 
of Eqs. (||) and (Q) in quite different ways. This prob- 
ably stems from the different role the block states \ip B ) 
play in matrix element storage, off diagonal with respect 
to total block occupation 0t s in singlc-fermion functions 
(if) B (^lf) \cj \ ip 3 ^ 3 )) but diagonal for the charge corre- 
lations (^ B (9lf) \n j \^ B ( ( R 1B )). One last comment con- 
cerns the q — charge correlations. We found no evi- 
dence for this contribution in our data possibly because 
of the specific quantity we chose to calculate in Eq. |[ 

We wish to point out an interesting observation we 
made on the raw (unaveraged) occupation rij near the 
ends. We can fit the occupation by the relation 



: + no cos(2k F j + 4>) (J) 



(8) 



where n = 0.25 and 0.25 < no < 0.4. A similar observa- 
tion has recently been reported in Ref. ^5[ The broken 
symmetry state resulting from the pinning at the chain 
ends forces the local occupancy variation (Srij), where 
Srij = rij — n, to be equal to the root mean square fluc- 
tuation y/ (SnoSrij). 

We wish to end this subsection by examining the sit- 
uation for Nb = 10. Why such a small number of block 
states? We have already stated that our calculations were 
made with Nb < 128. This is at the limit of our compu- 
tational capabilities. If the chains were independent, this 
would be equivalent to having some 10 block states per 
chain, which is not large indeed. At such small values of 
Nb , we had to introduce another coherence length £ c for 
the charge correlations, 

C2(r) = C cos(2/c F r + tp c ) 

x (£ c sinh(r/£ c )r 2i ^exp(-r/4) . (9) 



We used a fitting procedure which weighed more heavily 
the data for r < £ so as to be able to recover key pa- 
rameters with values close to those at Nb — 42. Trunca- 
tion errors run typically at the level of a few times 10~ 5 . 
This is considerably larger than for Nb = 42. Table || 
gives some of the parameters of the fit. We have used 
the values of A of Table I. a and K c compare favorably. 
The ground state excitation energy has appreciably in- 
creased and the coherence length has shortened. They 
are rather featureless, a signature of the small number 
of block states. Note that £ c ~ d c is somewhat "elastic" 
in the sense that its value can drift significantly without 
marked effect (within the 5 % error bar) on the fit. 



V. COUPLED CHAINS 

Now that we have acquired sufficient experience and 
confidence in data management, we can tackle the study 
of two coupled chains. One final observation is war- 
ranted. We found that convergence of our DMRG al- 
gorithm could only be achieved relatively quickly for an 
odd number of fermions Nf = (2iV/o+l)- This can be un- 
derstood in view of our discussion in Section [V A. First 



of all, fermions are pushed to the ends by the Coulomb 
repulsion. The fermion occupancies on both chains thus 
start site-synchronized at the ends. Secondly, the trans- 
verse hopping favors out of phase occupancies on the 
chains. Fig. ^ illustrates this. It shows the on-site charge 
(( n j,i) + ( n j,2)) and polarization — { n j.i)) that is 

typical of the broken-symmetry ground state. The larger 
V is, the more pronounced the out of phase character (po- 
larization) is and the larger the 4fc p charge component. 
This behavior can be achieved more easily in a state in 
which one chain has one extra fermion. The chain with 
(Nfo + 1) fermions is more compressed and cannot easily 
sustain solitons. But then, the chain with Nfo fermions 
can easily accommodate a kink-antikink pair which then 
allows the inner part of this chain to be out of phase with 
the ends and with the other chain. Broken symmetry 
states also have the advantage of focusing the computa- 
tional resources to a single non-degenerate state instead 
of splitting them between degenerate states thus decreas- 
ing the numerical coherence length. We have thus chosen 
to do our calculations for N — 150. Our truncation error 
varies from 5 x 10~ 5 for Nb — 42 to 5 x 10 -6 for the 
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FIG. 4: On-site charge (full line, + sign) and polarization 
(dotted line, - sign) in typical broken symmetry states at V = 
2, 6 and for t± = 0.5. 



larger values Nb = 128. We have kept the A values of 
the single chain. 

We shall first look at single-fermion behavior and then 
at some two-fermion correlation functions. 



A. Single-fermion transfer function 

The perturbative renormalization group formulation in 
Rcf . U presented a unified description of the renormaliza- 
tion of t± and the generation of interchain couplings in 
quasi-lD solids. Although the basic elements are present, 
the treatment is perturbative and subject to caution for 
large interactions. The bosonization approaches have the 
potential to do better in this respect since the single- 
chain interacting spinlcss fcrmion problem has an exact 
solution. 

Let us then examine the predictions of bosonization. 
There are two key treatments which look at our HamiL 
tonian from two different perspectives. Nersesyan et al.iB 
(NLK) bosonize the chaia-fermion operators Cj t p, whereas 
Yoshioka and Suzumuratil (YS) do so with the band op- 
erators 



J,2)/V2 



(10) 



Here a = ±1 is the band index. One has k± = (1 — 
<j)(n/2). The procedure yields two separated sectors, po- 
larization and occupation, that, by analogy to a chain of 
spins 1 in a magnetic field oc t±, are labelled spin and 
charge. 

In NLK, there is spin-charge separation and the 
Coulomb interaction is absorbed within the stiffnesses 
K c and K s . Here K c has the same value as for a single 
chain. tj_ appears in the spin sector and acts as the 
generator of interchain two-fermion couplings G and G 



corresponding to particle-hole and particle-particle pair 
hopping. The renormalization group (RG) equations 
for G reads G' = 2(1 - K S )G + (K s - K s )t 2 ,where 
t = A/(2nu s ), u s is the spin excitation velocity, 
K s = 1/K S , and the prime indicates the derivative with 
respect to £ = In (max (to / E ,T/ E ,vpk/ Eo)) where E 
is the starting energy scale. This is discussed in Refs. |l| 
and ^ The equation for G is obtained by the substitu- 
tion K s <=* K s . Note that G(£ = 0) = G(Qjj-e= 0. This 
RG equation is different from Yakovenko'st3 whose r 2 
term is larger by a factor 8ir 2 . This does not change the 
qualitative behavior of the equations, only the numbers. 
Furthermore, it is more in line with the coefficients of the 
RG equations in Ref. [| The RG also renormalizes t±, 
t' = (2 — A s )r where A s = \(K S + K s ), and K s is gov- 
erned by (In K s )' = \ (k 2 G 2 - K s G 2 ^j . Note that one 

has K s (£ = 0) = K c . There are additional contributions 
of order r 2 to this last equationu, a fact acknowledged 
by NLK, but which will remain unexplored by us. When 
G reaches strong coupling (G ~ 1) then a gap opens in 
the spin sector and only the charge sector contributes to 
power laws. 

In YS, there is also spin-charge separation for small t±. 
But the authors point out that this is no longer true for 
large transverse hopping, when the Fermi velocities for 
the two bands are appreciably different. This should be 
kept in mind as our values of t± > 0.1 should qualify. The 
charge sector behaves essentially as in NLK, with r\ = K c . 
The spin sector, however, transforms the Coulomb inter- 
action into two inter-band couplings g2<j>+ and g2<j>- such 
that the spin stiffness = 1 at I = 0. This dichotomy 
between spin and charge sectors is satisfactory for small 
V. The RG equations (when corrected for typographical 
errors) yield solutions qualitatively similar to NLK for 
t±(£). But it rapidly becomes annoying at large interac- 
tions since the spin sector quickly, too quickly, goes to 
strong coupling when (?2<^+ or <?2<£- ~ 1 and develops a 
gap. It then appears that the NLK approach should do 
better before the single-fermion dimensionality crossover 
and YS after the crossover, when bands have formed. 

Let us now focus on the predictions of these models 
for the transverse hopping as a function of length scale 
x(£) = exp(£ ). We first start with NLK at I = 0. The 
RG equations will switch to YS once the dimensionality 
crossover is reached, when t(£ x ) ~ (27r) _1 , via the map- 
ping proposed in NLK. We have chosen this value as it 
produces crossoj/ers at scales comparable with the ones 
of Bourbonnaisa and YS. With this choice, the renormal- 
ized transverse hopping 

t±(t) = (2nu s /A)T(£)exp(-£) « 2ttv f t(£)/x(£) 

leads to a dimensionality crossover at 

tx(4) = (27rti s /A)r(4)exp(-4) 

« u F exp(-4) = v F /x(£ x ) , (11) 

in which we have set A = 1 and u s = v F . This is the re- 
sult of Bourbonnais for the crossover. Fig. ^|(a) shows the 
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FIG. 5: Transverse hopping as a function of the length scale 
normalized to the bare value for various situations. The 
curves for V = 2 are for 0.5 > t± > 0.1 from top to bot- 
tom, (a) NLK bosonization results. The V = 6 curves are for 
all values of t±. (b) DMRG results. The V — 6 curves are for 
0.5 > tj_ > 0.2, from top to bottom, respectively. The stars 
are set at the values of the coherence length £. 



results for the renormalized transverse hopping t±(£) at 
both intermediate (V = 2) and strong (V = 6) Coulomb 
interactions. This was obtained using the K s = K values 
in Table | at I = 0. On the one hand, the V = 2 results 
are typical of small V and large t± and clearly show the 
dimensionality crossover. The "ladder" diagram is typi- 
cal of this crossover regime. The renormalized hopping 
after the crossover satisfies the relation 



(t±(4)/<±) « (*J 



W(l-a) 



(12) 



as proposed by Bourbonnais and verified by Capponi. 
On the other hand, the V = 6 curves for different t± all 
superimpose, a signature of the confined regime (a > 1). 

The DMRG calculations we now present were taken at 
Nb = 96. The following band transfer functions were 
calculated 



Cl(r,a) = C1((N -r)/2,r,a) 



(13) 



where C1(J, r, a) = 
Cl(r, c) was fitted 
be made by allowing 
bands a to vary as a 
single chains. This 
transverse hopping 
could then calculate 
from 



- (Cj +r , a Cj, a ). Each of these two 
with Eq. (0). The fit could only 
; the kp(r, a) which characterize the 
function of r, a situation not seen in 
is how the renormalization of the 
manifests itself most directly. We 
the renormalized transverse hopping 



tx{r) ~ [kp{r, cr = l) — fcf(r, a = —1)] /vf 



(14) 



The results for this last quantity are reproduced in Fig. 
|(b). The four upper curves for V = 2 clearly show the 
"ladder" characteristic of the crossover situation. This 
crossover, from Eq. (|ll|), is approximately at r x 
vp /t±(r x ) which has been reached before the coherence 



length is reached for tj_ > 0.2. The renormalized trans- 
verse hopping is also rather flat for r > r x . This is not the 
case for the t± = 0.1 curve which has a "drooping" (neg- 
ative slope) characteristic. This intermediate value of V 
has this nice property of showing both the crossover and 
its absence due to numerical coherence. Moreover, the 
predicted renormalization law of Eq. (|l^) gives a value 
a ~ 0.26 for the top four curves, a value in excellent 
agreement with the single chain one. This is somewhat 
accidental as the same analysis performed on the V = 1 
data for 0.2 < t± < 0.5 yields a value of a ~ 0.14, some- 
what larger that its value in a single chain. The problem, 
if there is one, may lie with Eq. (14) which does not ac- 



count for the change in Fermi velocity that occurs when 
the values of k± = 0, n are appreciably different (large 
t_i_) or due to renormalization. One must bear in mind 
that the DMRG is on a lattice and the fermion energy is 
far from the linear dispersion of the RG. Non-universal 
band edge effects (pre-RG) might be considerable and 
this might explain the discrepancy at the sizeable val- 
ues of t± we have used. We find no clear evidence for a 
crossover at V > 3 as £ < r x and the t± (r) curves do not 
flatten out at r x . 

The V = 6 curves are bunched together and show con- 
finement. This is in qualitative agreement with bosoniza- 
tion. Our results indicate that the coherence length sets 
the scale for the evolution of the transfer function. Fur- 
ther evolution, such as crossovers, are blocked at dis- 
tances beyond £. The splitting, albeit small, of the dif- 
ferent V = 6 curves is caused by just this effect. The 
coherence lengths are shorter for the larger values of t± 
thus leading to larger, that is less renormalized, values of 
t ± (£). ^ 

It is important to understand the connection between 
t_i_(r) and t±(£). t±(r) can be viewed, from Ref. ||, as 
an average over all momentum scales k < k r = lix jr. It 
is thus an average of t±(£) over £ < h\(v pk r / Eq) . Thus 
clearly t ± (r = 0) = t ± (£) f t ± (£ = 0). That is why 
t±(r = 0) is larger than t±(r) yet smaller than the bare 
value tj_. 

Fig. |^(a) shows the fitted values for a for both k± 
I ). 7r. They are surprisingly nearly independent of t± or of 
the existence or not of a dimensionality crossover. This 
is not expected from bosonization which would predict a 
variation of K s {£) and thus of a. This is perhaps a conse- 
quence of the absence of true spin-charge separation. YS 
also predict a change in exponent at the crossover. Fig. 
||(b) shows the average numerical coherence lengths for 
k±_ = 0, 7r. The coherence lengths monotonically decrease 
at constant V as t± increases. Curiously, £ is smallest at 
V = 2,3. We will come back to this in the next subsec- 
tion. 

Quite obviously, it would be futile to study the scaling 
behavior of these quantities with the chain length since 
£ is much shorter. But it would be possible to study the 
scaling with (Nb) 1 which acts like a temperature. In 
an attempt to better understand our results, we thus re- 
peated the calculations at other values of Nb- Fig. |?] 
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FIG. 6: Fitted values for the DMRG transfer function at dif- 
ferent transverse hopping and Coulomb interaction strengths, 
(a) Power law exponent results. The downward pointing tri- 
angles are for k± = n while the upward pointing ones are for 
k± — 0. (b) Mean values of the numerical coherence length. 



shows the scaling behavior of the ground state excita- 
tion energy and the reciprocal of the coherence length 
for V — 2. The gap is seen to extrapolate, as Nb — ► oo, 
to small values of E x compatible with N = 150. There is 
thus no indication of an energy gap. But no longer 
correlates well with E x in contrast to results for sin- 
gle chains. As a matter of fact, the extrapolated value 
increases quasi-linearly with t±. The same behavior is 
found for all values V > 1. Fig. || shows the Nb — > oo 
extrapolated inverse coherence lengths for V — 2 and 
V — 6 as well as the reciprocal of the threshold values to 
strong coupling coming from numerical solutions of the 
bosonization equations. We see there is a tight corre- 
lation between bosonization and DMRG behaviors with 
t±, though the length scales are different, in both sit- 
uations with or without a crossover. This suggests the 
coherence lengths are affected by the oncoming strong 
coupling regime at which our fitting formula, Eq. (|^), 
would no longer be valid. 



B. Two- particle correlations 

Following Eq. [|, we define the on-site charge [y = +) 
and polarization [y = — ) correlation functions C2(r, v) — 
C2{{N - r)/2,r, v), where 

C2(j,r,u) = ((n Vtj+r - (n Vj j+ r )) {n v>j - (n v>j ))) , (15) 

and 

n±,j = rij.i ± rij t 2 ■ (16) 

What should we expect for the power law exponents 

C2{r,v) cc r~ K » (17) 

characterizing these correlation functions? It is quite 
clear that as t± — > 0, the chains become independent and 
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FIG. 7: Ground state excitation energy E x (full symbols) and 
inverse coherence length £ -1 (empty symbols) as a function 
of Ng 1 for V = 2. The curves are for 0.5 > t± > 0.1 from 
bottom to top, for E x , and top to bottom for 
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FIG. 8: Inverse coherence length extrapolated to an infinite 
number of block states, circles, and reciprocal of the thresh- 
olds to strong coupling calculated from the bosonization equa- 
tions, squares, as a function of the transverse hopping. The 
results for V — 2, full symbols, and V = 6, empty symbols, 
are shown. 



one should have K v = K c + K s = 2K C . This is consistent 
with NLK. From the band perspective, the approach of 
YK predicts K + = (j] + rjj,) and K_ = (77 + 1/^0) when 
there is no gap. Since ry^ > 1 for repulsive interactions, 
one should expect K + > Jf_. 

We have calculated these correlation functions and 
tried fitting them to forms reminiscent of Eq. ([)]). The 
polarization correlation function behaves nicely and can 
be fitted with 

C2(r,-) = CLcos((A*(r,l)+*jr(r,-l))r + p_) 

x (£_ sinh(r/£-)r A '- exp(-r/cL) (18) 

using the filtering procedure explained in the data pro- 
cessing of single-chain correlations. The exponents K_ 
are essentially those of the charge correlations of single 
chains as witnessed in Fig. |9|. Again, this is unexpected. 
It is as if (77 + I/770) or (-Ke + Ks) remained nearly con- 
stant. This behavior is like the one observed for the sin- 
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FIG. 9: Power law exponents for charge, full circles, and po- 
larization, full boxes, for various values of V and t±. 



gle particle exponent a. This points to the breakdown of 
spin-charge separation for a highly non-linear dispersion. 
As for £_ and c?_ , we find they are generally larger than 
£, sometimes by a sizeable factor, and also fairly "elastic" 
as observed for single chains. 

The charge correlation function proved a bit more sub- 
tle to fit. The charge correlations generally decrease 
much faster than the polarization correlations. They 
rapidly become quite noisy when they reach the limits 
set by the truncation error. This can easily be moni- 
tored on a logarithmic plot of the amplitude of C2(r, +) 
as a function of r. There is a sudden break in the gen- 
eral linearlike decrease. The range of values in which a 
function having an amplitude of the type proposed in Eq. 
( |l8| ) can be appropriate is often limited to r < 30 when 
the coherence lengths are small. We have extended our 
calculations to Nb = 128 for a few cases as a check of 
the sturdiness of our calculations. It is mostly £_ that 
is affected by a noticeable increase. Curve fitting with 
many parameters can thus become delicate and leads to 
large uncertainties which we estimate at 10% for K + . We 
found that we could fit C2(r, +) with a form equivalent to 
Eq. ( |l8|) with v = + instead of v = — for all cases V > 3. 
The fast Fourier transform reveals only one broad u 2kp" 
wavenumber. K + is always of the same size or a bit larger 
than K— in this chain-like regime as can be seen in Fig. 
^. For those situations V < 2, the fast Fourier anal- 
ysis reveals two wavenumbers associated with 2kp(r, 1) 
and 2kp(r, — 1). We sometimes fitted the most prominent 
wavenumber or sometimes both, when the separation was 
not big enough, with a form that combines contributions 
of both bands (see Eq. |l^) 

C2(r,+) = ^2C +itT cos(2k F (r,(r)r + (p +<a ) 

<J 

x (£ + sinh(r/£+)r K+ exp(-r/d + ) .(l 9 ) 

As can bee seen from Fig. K + is appreciably larger 
than K_ in the low V regime. We found £ + and d + 
to be somewhat smaller than £ and still "elastic" . It is 



FIG. 10: Transverse hopping as a function of the length scale 
normalized to the bare value for V — 3 (full lines) and V — 4 
(dotted lines). The curves are for 0.5 > t± > 0.2 from top to 
bottom. 



quite apparent from Fig. ^ that there is a discontinuity 
in K + going from V = 2 to V = 3, with values drop- 
ping from 2 to 1. This behavior is, we believe, intrinsic 
to the coupled chains and is not a result of short coher- 
ence lengths since these are of the same size for V = 2, 3. 
It indicates a crossover from a single-fermion dominated 
behavior, when the charge fluctuations are propagated 
through independent band-like particle and hole motion 
showing two wavenumbers in the charge response func- 
tion, to a pair behavior, when the fluctuations are carried 
by particle-hole pairs having chain-like character shar- 
ing a common wascenumber. Such a crossover was de- 
duced by Capponll3 in the same range of values of V. 
In this crossover region, there is fierce competition be- 
tween single-fermion behavior and coherent pair motion. 
The computer resources are shared between the two com- 
peting behaviors and this explains the shorter coherence 
lengths in the crossover region. In spite of the domina- 
tion of the pair motion, the transverse hopping is not 
necessarily confined as Fig. [h] shows. There are no clear 
crossovers at the larger t±, using the criteria discussed 
in the previous subsection, but the hopping is consider- 
able and decreases only slowly with distance. There may 
well be coexistence between coherent interchain pairing 
and some sort of interband single particle behavior before 
the strong-coupling limit is reached, a possibility raised 
byNLK. 



It should finally be m entioned that, as was de- 
scribed in Section IV B for single chains, we find 



a broken symmetry state in which the total on- 
site occupancy variation (Sn + j) is proportional to 
the square root of the charge correlation amplitude 
(£+ smh(j /£ + ))~ K+ ^ 2 exp(— j/2d+) near the ends. 
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VI. SUMMARY 

Let us now sum up the more important findings on 
the quarter-filled double-chain problem. Our working hy- 
pothesis is that the fitting functions and the data treat- 
ment we have used for the analysis of one- and two- 
particle functions, which were validated for single chains, 
can safely be carried over to two coupled chains. 

1. The central result, which motivated this whole 
study, is the confirmation of strong renormalization of 
the inter-chain hopping in the presence of the Coulomb 
interaction V. This was measured by monitoring the dif- 
ference in the Fermi momentum of the two bands beyond 
the dimensionality crossover. The extent of the renormal- 
ization is much stronger that expected from the RG or 
bosonization at small V. It is considerable at the larger 
V. It is most likely going to zero as can be concluded 
by observing the trend in t±(r)/t± which decreases with 
t± . The proposed crossover law, Eq. ( |T^ ) , was seen to be 
only approximately valid when the crossover exists. It is 
probably our definition, Eq. (p^), and band edge effects 
that are responsible. 

2. The single- fermion transfer function power-law ex- 



ponent a and the polarization correlation function expo- 
nent K_ were seen to be essentially those of the single 
chain. This is interpreted as a sign of possible violation 
of spin-charge separation. 

3. The charge correlation function shows superposition 
of two wavenumbers and has an exponent that is 
twice as large as the one for polarization for V < 2. This 
indicates the existence of a crossover from independent 
particle-hole motion to correlated pair motion between 
V = 2 and V = 3. 

4. The ground state excitation energy, when extrapo- 
lated to Ng — > oo, is as expected for a finite system of 
150 sites. The coupled chains thus remain gapless. There 
are also indications, in the same limit, of possible strong 
coupling regimes as inferred from the linear-like t± de- 
pendence of the extrapolated inverse coherence lengths, 
much in the same way that bosonization predicts. 

The rather short numerical coherence lengths we have 
encountered with the DMRG, although expected from 
the study on single-chains, did put some stress on the 
curve fitting procedures. This resulted in fairly large er- 
rors in our exponents. 



* Electronic address: carou&pkysiquc.ushorb.ca] 
' Electronic address: cbourbon@physique.usherb.ca 

1 J. Solyom, Adv. Phys. 28, 201 (1979). 

2 Y. A. Firsov, Y. N. Prigodin, and C. Seidel, Phys. Rep. 
126, 245 (1985). 

3 C. Bourbonnais and L. G. Caron, Int. J. Mod. Phys. B 5, 
1033 (1991). 

4 C. Bourbonnais, in Les Houches, Session LVI (1991), 
Strongly interacting fermions and high-T c superconductiv- 
ity, edited by B. Doucot and J. Zinn- Justin (Elsevier Sci- 
ence, Amsterdam, 1995), p. 307. 

5 J. Voit, Rep. Prog. Phys. 58, 977 (1995). 

6 V. J. Emery, in Highly Conducting One-Dimensional 
Solids, edited by J. T. Devreese, R. E. Evrard, and V. E. 
van Doren (Plenum Press, New York, 1979), p. 247. 

7 S. Brazovskii and Y. Yakovenko, Sov. Phys. JETP 62, 1340 
(1985). 

8 H. J. Schulz, in Proceedings of Les Houches Summer School 
LXI, edited by E. Akkermans, G. Montambaux, J. Pichard, 
and J. Zinn- Justin (Elsevier, Amsterdam, 1995), p. 533. 

9 H. J. Schulz, G. Cuniberti, and P. Pieri, in Field Theories 
for Low-dimensional Condensed Matter Systems, edited by 
G. Morandi, P. Sodano, A. Tagliacozzo, and V. Tognetti 
(Springer, 2000). 

C. Bourbonnais and D. Jerome, in Advances in Synthetic 
Metals, Twenty Years of Progress in Science and Tech- 
nology, edited by P. Bernier, S. Lefrant, and G. Bidan 
(Elsevier, New York, 1999), pp. 206-261, arXiv:cond- 
mat/9903101. 



11 S. Biermann, A. Georges, A. Lichtenstein, and T. Gia- 
marchi, cond-mat/0107633. 

12 S. Capponi, D. Poilblanc, and E. Arrigoni, Phys. Rev. B 
57, 6360 (1998). 

13 M. Fabrizio, Phys. Rev. B 48, 15838 (1993). 

14 H. Yoshioka and Y. Suzumura, J. Phys. Soc. Japan 64, 
3811 (1995). 

A. A. Nersesyan, A. Luther, and F. V. Kusmartsev, Phys. 
Lett. A 176, 363 (1993). 

16 V. M. Yakovenko, JETP Lett. 56, 510 (1992). 

17 P. Donohue, M. Tsuchiizu, T. Giamarchi, and Y. Suzu- 
mura, Phys. Rev. B 63, 45121 (2001). 

18 S. R. White and R. M. Noack, Phys. Rev. Lett. 68, 3487 
(1992). 

19 S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 

20 S. R. White, Phys. Rev. B 48, 10345 (1993). 

21 I. Peschel, X. Wang, M. Kaulke, and K. Hallberg, 
eds., Density-Matrix Renormalization: a New Numerical 
Method in Physics, vol. 528 of Lecture Notes in Physics 
(Springer, Heidelberg, 1999). 

22 K. Yonemitsu, Synth. Met. 120, 845 (2001). 

23 F. D. M. Haldane, Phys. Rev. Lett. 45, 1358 (1980). 

24 H. Frahm and V. E. Korepin, Phys. Rev. B 42, 10553 
(1990). 

25 S. R. White, I. Affleck, and D. J. Scalapino, cond- 
mat/01 11320. 

26 The Pearson correlation coefficient is equal to the covari- 
ance of the two variables divided by the products of their 
standard deviations. 



